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ABSTRACT 

Internal (IS) or external (ES) shocks in relativistic expanding shells arc currently the 
best known mechanism for producing GRBs. We calculate the hydrodynamic condi- 
tions and the cooling processes in internal shocks, that occur when one layer overtakes 
another. We compare these conditions with the conditions in external shocks that occur 
when the shell encounters the ISM. We find that the complicated temporal structure 
observed in GRBs could be easily explained in IS scenario. Unlike ES, the observed 
temporal structure simply reflects the temporal behavior of the internal engine that 
produces the shell. We also find that opacity and in particular the opacity for pair 
production poses strong constraints on the parameter space and consequently bursts 
with very narrow peaks are expected to be optically thick and would not contain GeV 
photons. Finally we find that as in ES synchrotron is the most likely radiation process. 



Subject heading: gamma-rays: bursts-hydrodynamics-relativity 



1 INTRODUCTION 

A cosmological 7-ray burst (GRB) occurs, most likely, when 
the kinetic energy of a relativistic shell turns into internal 
energy and cools by non-thermal radiation process. The ki- 
netic energy of the shell can be transformed into internal en- 
ergy by decelerating onto the ISM (Meszaros & Rees 1992) , 
or by internal shocks in the shell (Rees & Meszaros 1994; 
Narayan, Paczyhski & Piran, 1992). 

Sari & Piran (1995, denoted hereafter SP) have dis- 
cussed the hydrodynamics of the deceleration of the rela- 
tivistic shell on the ISM and internal shocks. If the shell 
is thin and its Lorentz factor is low, internal shocks form 
before a considerable deceleration on the ISM takes place. 
Provided that the Lorentz factor of the shell is not uniform, 
with fluctuations of the order of a few, these internal shocks 
can extract a considerable fraction of the shell's kinetic en- 
ergy and produce the observed 7-rays. 

The internal shock (IS) scenario has several advantages 
over the external shock (ES) one. Mesazaros and Rees (1994) 
have pointed out that IS requires, quite generally, a lower 
Lorentz factor. This allows slightly larger amounts of bary- 
onic load in the fireball and it eases, somewhat, the con- 
straints on the source. We demonstrate in sections II and 
III that the hydrodynamics of IS also offers a natural expla- 
nation to the multiple time scales observed in the temporal 
profiles of the bursts. Such an explanation is missing in the 
ES scenario. 

The observed spectrum indicates that the source must 
be optically thin. Two processes are relevant for determining 
the optical depth: Compton scattering on the shell's elec- 
trons and pair production between the 7-ray photons. The 
relativistic motion of the shell releases the constraint that 



the optical depth poses to cosmological GRB sources. In fact 
this is the main reason to believe that highly relativistic 
motion is an essential ingredient of any cosmological GRB 
model (Goodman, 1986, Paczyhski 1986; Krolik and Pier, 
1989; Piran, 1996). For the ES scenario in which the kinetic 
energy is deposited at ~ 10 15 — 10 16 cm the optical thinness 
condition leads to a modest limitation on the Lorentz fac- 
tor 7 > 100 (Fenimore Epstein & Ho, 1993; Woods & Loeb 
1995; Piran 1996) which is not a serious restriction. For ES 
a stronger limit on 7 arises from the temporal structure. 
For internal shocks, in which the energy is deposited around 
~ 10 12 - 10 14 cm this poses an additional constraint which 
we examine in section IV. In the IS scenario the minimal 
Lorentz factor is determined by optical thinness condition. 
We show (in section IV) that the optical depth to pair pro- 
duction poses, usually a stronger constraint than the opti- 
cal depth to Compton scattering which Mesazaros and Rees 
(1994) used. 

Cooling must be rapid enough to allow for the shortest 
time scale. In a recent paper, Sari, Narayan and Piran (1996 
denoted hereafter SNP) have discussed cooling in the ES 
scenario. In sections V and VI we generalize this discussion 
to internal shocks, and we discuss the possible observational 
differences between cooling in the ES and the IS scenarios. 



2 TEMPORAL STRUCTURE 

Many bursts display a complicated temporal structure with 
numerous peaks (see e.g. Fishman and Meegan, 1995). We 
denote by tdur the total duration of the burst and by St 
the typical duration of the individual peaks. Few bursts are 
smooth with St m tdur - Most bursts are highly variable with 
St <C AT. Thus any GRB model must be capable of ex- 
plaining both time scales, the overall duration, tdur, and the 
duration of the individual peaks, St. We review here, briefly, 
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Figure 1. A photon emitted at (A) the outer side of a shell with 
thickness A will reach the oberve A/c before a photon emitted 
at the same time (in the obesrver's frame) from the inner side 
of the shell (B). A photon emitted at (A) will reach the observer 
R e /2^/ 2 c before a photon emitted by the same electron at 2r e (C). 
A photon emitted at (A) will reach the observer Re/"f 2 c before 
a photon emitted by matter that moves radially outwards at an 
angle I/7. 

several relativistic effects that determine the temporal struc- 
ture of the bursts (see Fenimore, Madras and Nayakshin, 
1996 for a detailed discussion). 

Consider a shell with a width A (in the observer frame) 
and a Lorentz factor 7 that converts its kinetic energy to 
internal energy and cools between R e and 2R e (see Fig. 1). 
A photon that is emitted at the outer edge of the shell will 
reach the observer A/c before a photon that is emitted at 
the same time (in the observer's rest frame) from the inner 
edge of the shell. Therefore: 

tdur > A/c . (1) 

We denote by 72 the Lorentz factor of the radiating 
matter (72 is always smaller than 7). A photon that is emit- 
ted by a given electron at R e will arrive R^j2^\c before a 
photon emitted by the same electron at 2R e (the difference 
arises from the the different time of flight of the photon and 
the electron from R e to 2R e ) (see Fig. 1). Finally the an- 
gular width of the observed emitting surface is I/72 due to 
relativistic beaming. A photon that it emitted radially out- 
wards will reach the observer R e J^j\c before a photon that 
is emitted by matter that moves at an angle I/72 from the 
line of sight (Katz, 1994) (see Fig. 1). The last two limits 
give rise to the same condition, which limits the duration 
of peaks. If the emission is continuous from R e to 2R e or if 
it spreads over a region with an angular width larger than 
Re/j2 then: 

St > R e /jic . (2) 

Within the ES scenario there are two possible regimes 
for deceleration of the shell on the ISM according to the 
value of the parameter £ (SP): 



£ = (Z/Af'VVs , (3) 

where / = (E/nismTripC 2 ) 1 ^ 3 is the Sedov length typically 
~ 10 18 cm. If £ > 1 the shock penetrating the shell is only 
mildly relativistic, the shell decelerates gradually, and 72 = 
7. We call this case the "Newtonian" scenario. If £ < 1 the 
shell is decelerated by one ultra relativistic shock, which 
reduces its Lorenz factor from 7 directly to 72 = 7£ 3 ^ 4 - This 
is the relativistic scenario. 

In either case the external shock satisfies: A < -Re/72 ■ 
Thus, if the energy conversion in the shell and the subse- 
quent emission is continuous or if the emission spreads over 
an angular size R e /y2 then tdur ~ St ~ dt. This corresponds 
to a smooth temporal profile with a single peak. Such pro- 
files are observed only in a small minority of the bursts. This 
conclusion is quite severe. It suggests that angular spreading 
will smooth all bursts and ES are not capable of producing 
bursts with complicated time profiles. 

A way out of the angular spreading problem, within the 
external shock scenario, is if the ISM is made out of small 
dense clouds with a large filling factor. There is an enhanced 
emission when the shell reaches ISM regions with higher den- 
sity. This requires large fluctuations in the ISM on the scale 
of » 10 12 cm (for peaks with St « 0.3 sec and 7 = 100). Sha- 
viv and Dar (1995,1996) suggest and explore in details the 
temporal structure for the scenario in which the shell inter- 
acts with a dense population of stars (in a globular cluster 
or a galactic center) which provide the desired high filling 
factor. Recently, Li & Fenimore (1996) have found that the 
distribution of peaks in long GRBs does not correspond to a 
random distribution. More specifically, the intervals between 
peaks has a log-normal distribution. Such a distribution is 
incompatible with a random distribution of objects that the 
shell encounters. We have already seen that the inner source 
cannot produce variability on such scale (because of the an- 
gular spreading). Hence, this feature, if confirmed, would 
poses a possible objection to the ES model. ^| 

The essence of the problem in the external shock sce- 
nario was the inequality: A < R e /72- We will see shortly 
that within the IS scenario we have, quite generally, A > 
-Re/72 an d 72 ~ 7- Consequently in this case A determines 
the overall duration, tdur, while the angular spreading as 
well as the radial size of the shell's fluctuations determines 
the duration of individual peaks, St. This, in turn, yields 
some interesting information about the inner source that 
drives GRBs. It requires that the source produces a shell of 
relativistic matter of width A in which fluctuations of or- 
der unity in the Lorentz factor should take place on a scale 
<5< A. 



* An alternative solution to the angular spreading problem would 
be if the source produces an irregular shell composed of many 
small blobs. Internal pressure will cause any small scale angular 
structure in the shell to spread radially and tangentially to size 
it/7 (Narayan and Piran, 1994, Piran 1994). However, recall that 
72 «C 7 in the relativistic ES scenario. If the shocked material 
cools before it expands from R e /-f to ite/72 it will produce a short 
peak on a scale t<jur/ c )72/7 (Sari and Piran, 1996). Temporal 
correlation in this substructure could be produced naturally by 
the internal engine 
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3 HYDRODYNAMIC OF INTERNAL SHOCKS 



Internal shocks occur when the relativistic ejecta from the 
inner source that drives a GRB is not moving uniformly. 
If some inner layer moves faster than an outer one it will 
take over. When the two layers collide two shocks appear 
- a forward shock propagating into the outer shell and a 
reverse shock propagating into the inner one. The forward 
and the reverse shocks have similar conditions behind them 
since it is most likely that both layers had (prior to the 
collision) similar conditions. This simplifies the discussion 
of internal shocks in which the two emitting regions have 
similar physical condition. In the case of external shocks we 
had to deal with each shock separately (SNP). 

We begin with a short discussion of the hydrodynamics 
of internal shocks in relativistic shells. Imagine a spherically 
expanding relativistic shell of thickness A in the observer 
frame moving with a high Lorentz factor 7. Assume further 
that the Lorentz factor varies by order of a few throughout 
the shell. Faster shells overtake slower ones. An inner shell 
will overtake an outer one at a radius, R s : 



(4) 



where S is the spatial scale of fluctuations in 7 within the 
shell, and we have ignored factors of order unity. A signif- 
icant fraction of the kinetic energy is converted to internal 
energy at this radius. 

The observed time over which internal shocks transform 
the relative kinetic energy of the two colliding layers (whose 
typical size is ~ 8) is of the order of 



Rah C = <V C 



(5) 



This corresponds to a single peak within a burst. Thus it 
can be identified with the observed quantity St. 

If more than one layer overtakes another one then other 
regions will collide and other shocks will form. These will be 
produced in different parts of the shell that are separated by 
radial distance of A. Therefore the total observed duration 
is given by 

tdur = A/c. (6) 
We define the ratio between S and A 

C^<1. (7) 

Comparison between Eqs. |, § and @ shows that £ is also the 
ratio between the duration of the individual peaks and the 
overall duration of the burst. 

The external shock which is produced by the ISM 
sweeps the shell by the time it arrives at radius 

Ra = / 3/4 A 1/4 , (8) 

(see SP). Clearly internal shocks are relevant only if they 
appear before the external ones. That is only if 



R s < Ra ■ 

Using Eqs. ^ and ^| condition ^| can be translated to: 

e /2 > c . 



(9) 



(10) 



This generalizes the result of SP who found that £ > 1 (for 
£ = 1) is a necessary condition for internal shocks. The con- 
dition can be turned into a condition that 7 is sufficiently 
small: 




Figure 2. Different scenarios in the A, (in cm) 7 plane for £ = 
5/ A = 0.01. Relativistic ES occur for large A and large 7 - upper 
right -above the § = 1 line (dark gray and light gray regions). 
Newtonian ES occur below £ = 1 - lower left - white region. IS 
occur, if there are sufficient variation in 7 below the £ = f 2 / 3 line 
(light gray and white regions). The equal duration t^ ur = lsec 
curve is shown for Newtonian ES (solid line) a relativistic ES 
(dotted line) and IS (dashed line). Note that a relativistic ES and 
an internal shock with the same parameters have the same overall 
duration t^ ur but different temporal substructure depending on 
5. 



7 < 280 C-V*( jkHHj-s/8, 



-3/8,3/8 



lOsec ' 



'18 



(11) 



We see that internal shocks take place in relatively "low" 7 
regime. Fig. 2 depicts the regimes in the physical parameter 
space (A, 7) in which different shocks are possible. It also 
depicts an example of a tdur =const. line. 

Provided that the different parts of the shell have com- 
parable Lorentz factor differing by factor of ~ 2, the shocks 
are mildly relativistic. The protons' thermal Lorentz fac- 
tor will be of order of unity, and the shocked regions will 
still move highly relativistically towards the observer with 
the approximately the initial Lorentz factor 7. Newtonian 
shocks have a limiting compression of factor of 7 (assuming 
an adiabatic index of relativistic gas, i.e., 4/3). The density 
behind the shocks is, therefore, comparable to the density 
in front of the shocks. In-front of the shocks the particle 
density of the shell is given by the total number of baryons 
E/jm p c 2 divided by the co-moving volume of the shell at 
the radius R s which is AtyR^A^. The density behind the 
shock is higher by a factor of 7. 

The hydrodynamical conditions behind the shocks are 
therefore: 



n = 7 E / Aivy c 4 rriptdurC 3 



(12) 
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4 OPACITY CONSTRAINS 



One of the advantages of the internal shock model is that 
it allows for short bursts with a lower value of the Lorentz 
factor 7. This eases the baryon purity constraints on source 
models. However, too low value of the Lorentz factor or too 
small emission radius could lead to a large optical depth. 
This poses another limitation on the parameters. 

The optical depth for Compton scattering of the pho- 
tons on the the shell's electrons at R s is given by: 



r e = {E/-y)/(4TvR 2 s )a T = 3.9 x 10" 5 C" 2 x 



(13) 



t-dur \ 

.wo) 



Er, 



V100/ 

The condition r e < 1 yields a lower limit on 7: 



7 > 13 



10 sec ) 



-2/5 

c 



-2/5 pi/ 



(14) 



which is not a strong limit. 

In addition, the radius of emission should be large 
enough so that the optical depth for 77 — > e + e~ will be 
less than unity (r 77 < 1). There are several ways to con- 
sider this constraint. The strongest constraint is obtained if 
one demands that the optical depth of an observed lOOMeV 
photon will be less than one (Fenimore, Epstein & Ho, 1993 
; Woods & Loeb). Following these calculations and using 
Eqs. ^ to express R s we find: 



7 > ioo( 



100 



-1/4 



(15) 



Eqs. 0, and |l| constrains 7 to a relatively narrow 
range. In fact the constraint due to 77 interaction is always 
more important than the constraint due to Compton scat- 
tering: that is r 77 > T e . In Fig. 3, we plot the allowed regions 
in the 7 and 5 parameter space. 

Three main conclusions emerge from the discussion so 
far. First, if the spectrum of the observed photons extends 
beyond lOOMeV (as was the case in the bursts detected by 
EGRET) and if those high energy photons are emitted in the 
same region as the low energy ones then the condition on 
pair production, r 77 , is much stronger than the condition on 
Compton scattering, r e . This increases the required Lorentz 
factors. In fact even if the spectrum has an upper cutoff of 
lOOMeV then it sets a lower limit of 7 > 100. Second, even 
the Compton scattering limit (which is independent of the 
observed high energy tail of the spectrum) poses a stronger 
than expected limit on 7 for bursts with very short peaks 
(see Fig. 3). The Lorentz factor should be as high as 10 3 to 
yield a peak of 0.3msec. Finally, one sees in Fig. 3 that only 
in a narrow region in the (S, 7) plane in which optically thin 
internal shocks are produced. The region is quite small if the 
stronger pair production limit holds and in particular there 
is no single value of 7 that can produce peaks over the whole 
range of observed durations. The allowed region is larger if 
we use the weaker limits on the opacity. But even with this 
limit there is no single value of 7 that produces peaks with all 
durations. The IS scenario suggests that bursts with narrow 
peaks should not have a very high energy tail and that very 
short bursts should have thermal spectrum. 




Figure 3. Allowed regions for internal shocks in the 5 (in cm), 7 
plane. Note that the horizontal 8 axis also corresponds to St, the 
peak duration by a devision by c. Internal shocks are impossible 
in the upper right (light gray) region. The lower boundary of this 
region depend on £ = 5/ A are marked by two solid curves, the 
lower one for £ = 1 and the upper one for £ = 0.01. Also shown are 
t 77 = 1 for an observed spectrum with no upper bound (dotted 
line), r 77 = 1 for an observed spectrum with an upper bound of 
lOOMeV (dashed line) and r e = 1 (dashed-dotted). The optically 
thin internal shock region is above the r = 1 curves and below 
the § = C 2 / 3 (solid) lines. 



5 SYNCHROTRON COOLING 

Following SNP we use the physical conditions derived in the 
previous section to estimate the cooling of the shocked re- 
gions. In order to estimate the radiation processes we need 
the energy density, e, the particle density, n, given by equa- 
tion [l^ as well as the magnetic field strength, B, and the dis- 
tribution of electron Lorentz factors j e . These last two quan- 
tities are difficult to estimate from first principles. Both de- 
pend on the microphysics of particle acceleration and mag- 
netohydrodynamics within the shocks. Our approach, fol- 
lowing SNP will be to define two parameters, e_e and e e to 
incorporate these uncertainties. We then constrain the val- 
ues of these parameters by requiring the model predictions 
to resemble the observed features of GRBs. 

The dimensionless parameter en measures the ratio of 
the magnetic field energy density to the total thermal energy 
density, e: 



es = 



U B 



87re 



(16) 



Using Eqs. and |l^ we can express B in terms of es and 
other observed quantities: 
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B = (720 G^C 1 (^) Eil\^r 3 . (17) 

The second parameter e e measures the fraction of the 
total thermal energy e which is in the form of the electrons' 
random motions: 

Ue 



Since the electrons obtain their random motions via shock- 
heating, we make the standard assumption that they develop 
a power law distribution of Lorentz factors, 



iV(7e) ~ if for 7 e > 7 e 



(19) 



Following SNP we fix the index f3 » —2.5 and we find 
that the average electron energy density is 



(20) 



Using Eqs. [l^ and [20] one finds that the electrons behind the 
shocks satisfy 



7e 



610e e 



(21) 



The value of 7 e ,min is sufficient for the calculations of the 
cooling time scale and efficiency that follow. The exact de- 
tails of the 7e distribution are not needed. 

We can estimate now the Lorentz factor of an elec- 
tron that is radiating an energy hv « lOOKeV: 



m e chv bs 
hv„ hs \ 1/2 



1/2 



= 1.1 x 10 V 



(tdur \ 
Ws) 



3/4 



e: 



; l/4 c l/2 > 
1/4 / 



(22) 



ViooKevy Uos/ Jl uoo/ 

We first check that electrons with this Lorentz factor are 
available. This requires, of course, fmin < which corre- 
sponds to 

!/2 ft, x3/4 

Viooy 



e e < 18ei 



1/41/2/ hu obs 1/4/ 7 \ ,„„, 

4 VlOOKevJ VlOsJ 



These electrons are always available. A similar situation oc- 
curs in a reverse external shock. But in a forward external 
shock for short burst durations (which force high Lorentz 
factor) the minimal electron energy is too high and most of 
the radiation emerges in higher energies than the observed 
band. 

Using the value of -y e we can estimate the synchrotron 
cooling time. The electron's energy divided by the power 
of its synchrotron radiation (with the appropriate Lorentz 
transformation from the electron's frame to the observer's 
frame) : 

-1/2 



ts 



100 KeV, 

4 /+ . \ 9/4 



(24) 



100/ VlOs/ 5 



3/4 



This cooling time is shorter than the cooling time es- 
timated for external shocks (SNP) since the magnetic field 
here is higher. It sets a lower limit to the variability time 
scale. Clearly, the burst cannot possibly contain spikes that 
are shorter than its cooling time. The cooling time, t syn 
is proportional to {hv) -1 / 2 as in the case of ES and in 
agreement with observations. This is a general feature of 
synchrotron radiation. An observation of this feature can, 



therefore, verify that synchrotron radiation takes place in 
GRBs but it would not distinguish between the ES and the 
IS scenarios. 

Using Eqs. ^ we calculate the duty-cycle D defined as 
the ratio of the cooling time to the total duration: 

'1/2 



(18) D = — ^ 

^dur 



tayn = 1.4 x 10~ 4 e 



/ hu obs \ " 
U00 KeV/ 



-3/4^3/2 

^ MOO KeV 



(25) 



3/4 



Uooy U0s, 

The duty cycle increases with burst duration, but it is short 
enough to explain the high variability even for bursts as 
long as lOOsec . This is different from the ES scenario in 
which the duty-cycle was only weakly decreasing with burst 
duration. Since duration changes by about four orders of 
magnitude it might be possible to distinguish between the 
two models. Note, however, that in internal shocks the syn- 
chrotron cooling is so rapid that another process (such as the 
acceleration within the shock) could easily be slower and it 
would determine the shortest observe time scale. Thus, the 
strong dependence of t 3y n on tdur might be screened and 
undetectable. 

The cooling time scale increases as the magnetic field 
decreases. Most GRBs are highly variable with a duty-cycle 
of order of 5% or less. This sets a lower limit to the value of 



^ o o „ / hv bs 

£B * 3 - 8 X 10 C llOOKeV 



-2/3 



(26) 



V 0.05/ 



D , 7 ye/s (U , 5/3 

llooJ \wl) ■ 

Unlike ES, the magnetic field here can be quite far from 
equipartition. This eases somewhat the requirement on the 
processes that build up the magnetic fields in the shocks and 
it is possible that the magnetic field dragged along the fluid 
from the source might be sufficient. 



6 INVERSE COMPTON SCATTERING 

Inverse Compton (IC) scattering may modify our analysis 
in several ways. IC can influence the spectrum even if the 
system is optically thin (as it must be) to Compton scat- 
tering (see e.g. Rybicki & Lightman, 1979). The effect of IC 
depends on the Comptonization parameter Y — 7 2 r e . It can 
be shown (SNP) that Y satisfies: 



Y = e e /e B 



C-B 



if e e < e B 
if e e > eb- 



(27) 



IC is unimportant if Y < 1 and in this cases it can be 
ignored. 

If Y > 1, which corresponds to e e > eg and to Y = 
v/ £e/es then a large fraction of the low energy synchrotron 
radiation will be up scattered by IC and a large fraction of 
the energy will be emitted via the IC processes. If those IC 
up scattered photons will be in the observed energy band 
then the observed radiation will be IC and not synchrotron 
photons. If those IC photons will be too energetic, then IC 
will not influence the observed spectra but as it will take a 
significant fraction of the energy of the cooling electrons it 
will influence the observations in two ways: it will shorten 
the cooling time (the emitting electrons will be cooled by 



© 0000 RAS, MNRAS 000, 000-000 



6 Sari & Piran 



both synchrotron and IC process). Second, it will influence 
the overall energy budget and reduce the efficiency of the 
production of the observed radiation. We turn now to each 
of this cases. 

Consider, first, the situation in which Y > 1 and the IC 
photons are in the observed range so that some of the ob- 
served radiation may be due to IC rather than synchrotron 
emission. This is an interesting possibility since one might 
expect that the IC process to ease the requirement of rather 
large magnetic fields that is imposed by the synchrotron pro- 
cess. We show here that, somewhat surprisingly, this is not 
the case. 

An IC scattering boosts the energy of the photon by a 
factor 7e. Therefore, the Lorentz factor of electrons radiat- 
ing synchrotron photons which are IC scattered on electrons 
with the same Lorentz factor and have energy hv in the ob- 
served range is the square root of % given by Eq. These 
electrons are cooled both by synchrotron and by IC. The 
latter is more efficient and the cooling is enhanced by the 
Compton parameter Y. We can now estimate the duty-cycle 
(the ratio of the cooling time to the overall duration): 

D = 0.0146" 3/8 6 e - 1/2 C 7/4 (^) x (28) 

/ tjur \ 13/8 E -r/8 ( hu obB X- 1 / 4 
UOsec/ 51 UOOKeV/ 

This duty-cycle depends strongly on the duration, a fact 
that was not observed. Additionally, it scales with the ob- 
served photon's energy as (fei/)~ 1//4 while observation are 
more compatible with (hu)^ 1 ^ 2 (Fenimore et. at, 1995). It 
is unlikely therefore, that the observed radiation is due to 
IC process. In addition, inspection of Eq. ^ shows that the 
magnetic field cannot be very small even in this case. If en 
is too small D will become larger than unity. 

Finally, if Y > 1 IC will influence the process even if 
the observed photons are not produced by IC. It will speed 
up the cooling of the emitting regions and shorten the cool- 
ing time, t S yn estimated earlier by a factor of Y. With this 
increase in the cooling rate the duty-cycle limit on the mag- 
netic field become extremely small. However, these very low 
magnetic field values are impossible since in addition IC also 
reduces the efficiency by the same factor, and the efficiency 
becomes extremely low as described below. 

The efficiency of a burst depends on three factors: first 
only the electrons' energy is available. This yi elds a factor e e - 
Second, there is an additional factor of of yj eg /e e if the IC 
radiation is not observed. Third, only the energy radiated by 
electrons with 7 > 7 is observed, and assuming a a power low 
electron distribution with an index (3 = —2.5 (see SNP and 
the previous section) this gives a factor of ("fmin/'y) 1 ^ 2 - The 
total efficiency is the multiplication of those three factors 
and assuming that es < e e it is given by: 

e tot > OMeTe^Ell* ( Jl) ^ ^ . (29) 

The efficiency can be rather high provided that the electrons 
and the magnetic field energy density are close to equipar- 
tition. This efficiency equation sets a stronger lower limit 
for the magnetic field than the duty-cycle lower. As in the 
ES case the efficiency increases when the burst duration de- 
creases. 



7 DISCUSSION 

We have calculated the hydrodynamic and cooling time 
scales of GRBs produced by internal shocks. We have found 
that internal shocks could easily explain the observed com- 
plicated temporal structure of GRBs: it directly reflects the 
temporal behavior of the inner engine that drives the GRB. 
Recall, that this is not the case in the ES scenario (SP). The 
duration is dictated by A, the width of the relativistic shell, 
which in turn corresponds to the total duration of the emis- 
sion of the internal source. The duration of individual peaks 
is determined by S the length scale over which the condi- 
tions in the shell vary significantly. The length 8, is again 
dictated by the temporal scale of variability of the internal 
source. This scenario agrees, therefore, with the conclusion 
of Li & Fenimore (1996) who discovered that the mid-peak 
intervals has a Log-Normal distribution which indicates a 
causal relation between peaks and concluded that the tem- 
poral structure is driven by the inner source. 

The opacity, and in particular the opacity to pair- 
production poses a sever constraint on the parameter space. 
The strongest constrain appears if we demand that an ob- 
served lOOMeV photon can escape freely from the system 
(Fenimore, Epstein and Ho, 1993; Woods and Loeb, 1995). 
The fact that there are such observed photons and that 
EGRET observes even higher energy photons suggests that 
this demand is reasonable. If the observed spectrum is un- 
bounded than this yields a strong limit and only a narrow 
strip is allowed for optically thin IS in the 7, 8 plane. The 
most worrisome fact is that there is no single value of 7 (or 
8) which can account for all kinds of observed bursts and 
the shape of this allowed region demands a correlation (in 
fact anti-correlation) between 7 and 8. This is not impossible 
but it seems unlikely. Alternatively this suggests that bursts 
with short peaks will be optically thick to lOOMeV photons 
and such photons won't be observed in those bursts. 

If there is an upper cutoff E max to the observed spec- 
trum then the optically thickness problem is weaker and 
pair-production simply set a limit of a minimal Lorentz fac- 
tor of ~ E max /2m e c 2 which is 7 ~ 100 for lOOMeV photons. 
Even in this case the optical depth to Compton scattering 
becomes larger than unity for small values of 8 - which cor- 
respond to very short peaks. The allowed region is wider but 
once more either there is some correlation between 7 and 8 
or we should observe optically thick narrow peak bursts. 

Like in the ES scenario the radiation process is most 
likely synchrotron. In spite of the higher energy density in 
the internal shock region (higher compared to the densities 
in the ES scenario) the magnetic field values required are, 
close to the equipartition value. Otherwise, the efficiency 
will become too low! Inverse Compton (IC) is possible - but 
quite unlikely. Somewhat surprisingly IC cooling requires 
magnetic fields that are as high as those required for syn- 
chrotron cooling. Again if these values are too small then 
efficiency constraints will rule out the model. 

Synchrotron cooling in the IC scenario suggests a cor- 
relation between the cooling time and the overall duration: 
t sy n oc i^ur- This is a high power and if t syn determines the 
duration of the shortest peaks it might be observed since tdur 
varies by four orders of magnitude. This should be compared 
to the cooling time in the ES scenario which varies like t 3 ^ 4 . 
However, the cooling rate in the IS scenario is so rapid that 
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this possible correlation may be masked by another process, 
e.g. the acceleration within the shocks, that would be slower 
and would determine the duration of the shortest peaks. 
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